Preserving Quantum States Using Inverting Pulses : A Super-Zeno Effect 
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We construct an algorithm for suppressing the transitions of a quantum mechanical system, ini- 
tially prepared in a subspace V of the full Hilbert space of the system, to outside this subspace by 
subjecting it to a sequence of unequally spaced short-duration pulses. Each pulse multiplies the 
amplitude of the vectors in the subspace by —1. The number of pulses required by the algorithm to 
limit the leakage probability to e in time T increases as T exp[\/log(T 2 /e)], compared to T 2 e _1 in 
the standard quantum Zeno effect. 
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Quantum computation is known to be more powerful 
than classical computation. Two well-known examples 
where quantum computation gives a clear advantage over 
the classical are searching and factorization. A key dif- 
ference is that, only in the quantum case, one can use the 
possibility of destructive interference of quantum ampli- 
tudes to suitably design algorithms where one can cancel 
out the amplitudes of undesired states. We want to apply 
this general strategy to the problem of state preservation, 
i.e. how to limit the system's evolution to a desired sub- 
space [lj. 

Several different approaches to this question have been 
studied. Error-avoiding codes depend on existence of 
subspaces free of decoherence due to special symmetry 
properties, and error-correcting codes Q depend on mon- 
itoring the system and conditional feedback to achieve 
this. Another important strategy for state preservation 
is to use the well known quantum Zeno effect, stud- 
ied by von Neumann, and others viz. preserving a 
given quantum state by very frequent measurements. In 
fact, the reduction of wavefuction is not essential for the 
physics of the quantum Zeno effect, and one can achieve a 
similar suppression of transitions by suitably coupling the 
system strongly to an external system for short intervals 
in the bang-bang control p| and dynamical decoupling 
strategies |6j , where the time evolution is always unitary, 
without any reduction of the wave-packet 0- 

In this paper, we develop a bang-bang control type 
algorithm to preserve quantum states that involves 
unitary kicks interspersed with evolution according to 
the system-environment Hamiltonian for unequal time- 
intervals between between kicks. The unitary kicks we 
employ are inverting pulses similar to those used in 
Grover 's quantum search algorithm The proposed 
algorithm does not require any symmetry properties and 
is much more efficient than the quantum Zeno effect. We 
will show that in our algorithm, the number of pulses 
required to keep the the quantum system in the same 
subspace up to time T with probability greater than 
1 — e increases only as T exp[^/log(T 2 /e)] , for large T, 
and small e, whereas it varies as T 2 e _1 and T 2 e~ 1//2 for 
the quantum zeno effect with and without measurement. 



We shall call the preservation of state using such non- 
periodic pulse sequences the super-Zeno effect. 

We consider a quantum mechanical system described 
by a finite dimensional Hilbert space TL, which is a direct 
sum of two orthogonal subspaces V and Q. We assume 
that the system Hamiltonian H is bounded. The unitary 
operator corresponding to evolution for a time interval 
t is given by Uo(i) = exp[— iHt]. When the system is 
subjected to a very short duration external fields pulse, 
we assume that the effect of the pulse can be represented 
by a unitary operator J. If the system is initially in the 
state after being subjected to the pulse, its state is 
J|</>>. 

A sequence of N pulses is specified by (N + 1) real 
numbers {xj},j — 1 to N + 1, where Xjt is the time- 
interval between the jth and (j + l)th pulse for j = 1 to 
N — 1, and XN+it is the interval between the last pulse, 
and the measurement of the state of the system. Here 
t is time interval between the initial preparation, and 
final measurement, and Xj = 1. In this case, the 

evolution operator is 



Wjv(t) = U (x A r +1 <)J . . . 3U (x 2 t)3V (x 1 t). (1) 

Our aim is to choose a sequence of the time-intervals 
between pulses given by {xj}, such that if the system is 
initially in a state \ip) G V , then the probability that the 
system is found in a state in Q after the pulse sequence 
is minimized. 

In this paper, we shall discuss only the case where 
J is an inverting pulse: J = Q P, where P and Q 
are the projection operators for the subspaces V and 
Q respectively. Clearly, J\tp) = —\ip), if \ip) G V, and 
3\ij}) = +\ip), if \ip) G Q. We first give an explicit 
construction of a recursively defined sequence U m (i) 
(see Eq. (4) below), such that the transition amplitude 
(q\U m (t)\p) is 0(t m+1 ). We then obtain quantitative 
bounds on the leakage probability, and compare the per- 
formance of our algorithm with the standard quantum 
Zeno effect. 

The basic idea of using inverting pulses to produce de- 
structive interference of quantum mechanical amplitudes 
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and reduce the transition rate from a particular subspace 
to others is quite straight-forward. Let \p) be a general 
state in V, and \q) a general state in Q. Let U(i) be a 
unitary operator satisfying U(t) = I + O(i) for t — > 0, 
where I is the identity operator and the matrix elements 
U qp , and U P q are of 0(t r ), with r is a positive integer. 
Then in the transition amplitude (g|U(t)JU(t)|p) a pre- 
cise cancellation of the amplitude of 0(t r ) term occurs, 
and one gets 

(q\U(t)JXJ(t)\p)=0(t r+1 ). (2) 

Similarly, it is easy to see that if V(t) is an operator with 
V(t) = J + 0(f), then again a destructive interference of 
amplitudes to leading order occurs in (<7|V(i) 2 |p), and if 
V(t) qp = 0{t r ), then, 

< (Z |V(£) 2 b)=0(f+ 1 ). (3) 

Using this result, it is straight forward to construct a 
pulse sequence where the transition amplitude is 0(t r ) 
for any positive integer r by recursion. We define opera- 
tors U m by the recursion relations 

U m+1 (f) = U TO (t/2) J U TO (t/2), for to even, 

= U TO (t/2) U m (i/2), for m odd. (4) 

with U = U (i). Then, by induction, it follows that the 
transition amplitude (q\U m \p) is of order 0{t m+1 ). 

We can write down U m (f) explicitly as a product of 
U (f/2 m )'s and J's. For example 

Ui(i) = U (f/2)JUo(t/2), 
U 2 (i) = [U (f/4)JU (f/4)] 2 , 
U 3 (i) - [U (f/8)JU (f/8)] 2 J[U (f/8)JU (f/8)] 2 . (5) 

Let jV m be the number of pulses used the sequence 
U m , then it is easily verified that 

N m = (2 m+1 - 2)/3, for to even, 

= (2 m+1 - l)/3, for m odd. (6) 

We can derive quantitative bounds to estimate the 
probability of leakage out of the subspace V for this se- 
quence U m . We will show that the leakage probability 
L m from subspace V to subspace Q for operator U m ap- 
plied for a total time T is 

L m = J2 l(?|U m |p)| 2 < 2 - m ^[ET] 2m + 2 . (7) 

Here -E = |H|, and we have used the definition of the 
norm of an operator A, given by 

\A\= S J \\A\M/\\\M (8) 

By making m large, L m can be made as small as we 
please. This becomes a state preservation algorithm 
when the subspace V consists of just the initial state. 
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FIG. 1: The plot of the upper bounds to the leakage proba- 
bility versus the number of pulses. The curves are marked Z 
for the Zeno with periodic measurement, U for unitary evolu- 
tion with periodic pulses, and S for the Super-Zeno cases, for 
ET — 1. Here T is the total time for evolution, and E is the 
norm of the hamiltonian. 



Note that the minimum time interval between pulses de- 
creases as 2~ m T 7 and this may limit the maximum value 
of to that can be realized in an experimental setup. 

For the TO-th sequence, , the leakage probability varies 
as 2~ m , while the number of pulses increases as 2 m . 
Thus, the leakage probability varies as exp[— (log N) 2 ] . In 
Fig. 1, we have plotted the upper bound on leakage prob- 
ability L m as a function of the number of pulses N m . For 
comparison, we have also shown the corresponding bound 
on the leakage probability for same number of pulses for 
the quantum Zeno effect, both with repeated measure- 
ments, and with unitary evolution with periodic inver- 
sion pulses. For the last two, these bounds are (ET) 2 /N 
and (ET^/N 2 rspectively ( see discussion below). The 
improvement in the super-Zeno over the others is quite 
large, even for fairly small values of to; e.g. for ET = 1 
and to = 4, which corresponds to only 10 pulses, the 
leakage probability is already of order 10~ 4 . 

To prove the inequality J7J), define 

Mnr 

Bmit)= \p)\q) K^WWIHQU^PI (9) 

We have already seen that B m is 0{t m+1 ). Since the 
only available timescale is 1/|H|, and B m is bounded by 
1 for large t, we can choose constant K m such that 

B m (t) < K m (Et) m+1 , for all t. (10) 

The proof of Eq. (J7J) is complete if we can show that 

K m < 2-™( m+1 )/ 2 . (11) 

This can be proved using the following simple mathemat- 
ical inequalities valid for all mutually orthogonal projec- 
tion operators P and Q, and all unitary operators U and 
V, with J = Q P (proof omitted): 

|QUJUP| < 2 |QUP| Min[l, |U-I|], (12) 
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and 

|QVVP| < 2 |QVP| Min[l, |V- J| ]. (13) 

Putting U = U r _i, with r odd, and V = U r _i with r 
even, then using |U — I| < Et and |V — J| < Et in the 
inequalities above, we get 

K r < 2- r K r _ 1 . (14) 

This, together with Kq — 1, implies Ea. i|ll|) . and com- 
pletes the proof of Eq.Q. 

We now compare the performance of our algorithm 
with that of state preservation using the standard Zeno 
effect. Let C(T,e) be the minimum number of pulses 
required to keep a quatum system in a precribed state 
for time T, with the leakage probability less than e. It 
seems reasonable to take this as the measure of the cost 
of algorithm. 

Recall that in the standard quantum Zeno effect, start- 
ing with an initial state \p), and measuring the projec- 
tor \p) (p| repeatedly, at times T/N, 2T/N, ■■■ , T, the 
probability of finding the system in the state \p) in each 
of these measurements is S(T) = \{p\e~ lIiT ' N \p)\ 2N . For 
small T/N, we can approximate this as [1- A(T/N) 2 ) N < 
(l-AT 2 /N), where A = (p\H 2 \p) - (p\H\p) 2 < E 2 . Thus 
the leakage probability is has an upper bound (ET) 2 /N , 
and C(T, e) varies as T 2 /e. 

If we use N equispaced unitary kicks J in time T, the 
cost C(T, e) varies as T 2 /i/e. This follows from the fact 
that the magnitude of the transition amplitude from sub- 
space V to Q after two pulses is < (ET) 2 /N 2 ), and , the 
magnitude of the transition amplitude after N pulses, 
and not the transition probability, has an upper bound 
(ET) 2 /N. 

In our case, the number of pulses N m varies as 2™, and 
the leakage probability varies as exp(— to 2 ), for large to. 
Hence, for fixed T and large to, the leakage probability 
L m decreases as 2~( log27V ' . Equivalently, for fixed T,a s 
e is decreased, C(T, e) increases as exp(ATy / log(l/e)), 
where K is a constant. This is slower than any power 
law in e. 

Now we discuss the behavior of C(T, e) for large T. In 
this case, Eq.JTJ) implies that C(T, e) increases at most as 
T 2 for large T. However, for large T, the bound Eq.JJJ) 
is too weak, and we can prove the stronger result 

C(T, e) < KET 2^°^^, (15) 

where K is some constant. 

We indicate the proof here. Consider 2 r+1 > ET > 2 r , 
where r is a non-negative integer. Assume to > r. From 
the first halves of inequalities (|12J) and (|13J) . B m (T) < 
2 r B m - r (T/2 r ). Then using the bound for B m - r , we 
get 

B m (T) < 2 r (ET2- m ) m - r+1 2 tm ~ rH *~ r+1} . (16) 

On taking logs, the inequality B m < y^e, becomes a 
quadratic inequality for to, which is easily solved to give 



to > log 2 (£T) + Vlog 2 (£ 2 T 2 /e) + 0(1), (17) 

for ET ^S> 1. Since the number of pulses increases as 2 m , 
we get Eq.JDjJl. 

Our result is similar to that found by Khodjasteh and 
Lidar 9] , who also used a recursive construction of pulse 
sequences, and found it to be better than the periodic 
pulse sequence. They considered pulse sequences for 
preservation of quantum states which work for all ini- 
tial states, and the pulse sequence is independent of the 
state to be preserved. However, they require more than 
one type of pulses to achieve this. In contrast, we allow 
only one type of pulses, and can only minimize the leak- 
age probability out of a subspace, and the pulse depends 
on the subspace to be preserved. 

We note that in our construction, the intervals between 
pulses take only two values, one of which is twice the 
other. This may be convenient in actual implementa- 
tion, as one can start with a source that generates peri- 
odic field pulses of the desired type, and then selectively 
block some of the pulses. There is no problem in relax- 
ing the assumption about pulse- width being infinitesimal. 
In fact, if the inverting pulse is error-free, then one can 
just use a periodic repetition of inverting pulses to get a 
perfect preservation of the state. Effect of noise in the 
inverting pulse would strongly affect the performance of 
the algorithm, and remains to be studied. 

It is possible to decrease the required number of pulses 
significantly if we allow the intervals between pulses to 
be varied continuously independent of each other. In 
general, one can try to find non- negative (N + 1) real 
numbers {xj}, such that for the matrix Wjv(t) defined 
in Eq.(l), the transition amplitude (q\Wisf(t)\p) is of a 
specified order 0(t m+1 ), as t — > 0. This is certainly pos- 
sible, by our explicit construction, for N = N m . Also, if 
this is possible to do for one value of N = r, it is also 
possible for N = r + 1, as one can always add the (r + 1)- 
th pulse at the end of time T, with x r +2 = 0. Hence 
we have proved that for any to > 0, there is an integer 
N m i n (m), such that for all N > N m i n (m), one can find 
an N- pulse sequence Wat(£) such that the transition 
amplitude (q\W N (t)\p) is 0(t m+1 ), for small t. 

Clearly, N m i n (m) < N m . We can improve this 
bound using a variation of the recursive construction 
of reflection-symmetric sequences ( i.e. those for which 
xj = xjy + 2-j, for all j) due to Yoshida 10]. For such 
sequences, we have WAr(t)Wjv(— t) = I for all t. This 
implies that if (<?|Wjv(t)|p) = 0(t r ), then r must be odd. 
Then construct 

W'(t) = W N (at)3W N (f3t)3W N (at), (18) 

with a = 1/(2+2"), j3 = 1 — 2a. Then it is easily checked 
that the 0(t r ) term in (q[W'\p) vanishes by construction. 
Hence the leading term in (<7|W|p) must be the next odd 
term, i.e. (D(t r+2 ). Starting with number of pulses N r = 
for r = 1 , we can recursively construct a sequence where 
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the transition amplitude is 0(t 2m+1 ) uses only 3™ 1 — 1 
pulses. Thus, we get N min (2m) < 3 m — 1, which is a 
significant improvement over N m i n (m) < N m . 

Further improvements are posible. For a given N, one 
can express (q|Wjv(i)|p) as Taylor series in powers of 
t. The coefficients of different powers of t are sums of 
matrix-elements of the type (q\H. ni 3H n2 3..\p), with co- 
efficients that are polynomials of {xj}. We try to set all 
such coefficients for powers of t up to to equal to zero, and 
solve the resulting polynomial equations for {xj}. This 
calculation is straight forward for small to. For m = 0, 1 
and 2, the sequences we get are Uo,Ui and U2. For 
larger to, U m or the W arc not optimal. 

The number of polynomial equation for equating to 
zero all terms up to order to are (2 m — 1). Interestingly, 
one can find a much small number of variables {xj} that 
do satisfy all of these. In particular, for to = 5, we can 
find six positive real numbers Xi, . . . Xe that satisfy all the 
31 equations generated. For still higher to, the equations 
become rather messy to analyse, even with Mathematica. 
However, in the cases that we could solve, the final solu- 
tion is surprizing simple. We omit the details, only men- 
tion here that for m = 4, the optimizing pulse sequence 
is {xj} = {/?, 1/4, 1/2 -2/3,1/4,/?} with j3 = (3->/5)/8. 
For m = 5, it is {7, 1/4 - 7, 1/4, 1/4, 1/4 - 7,7}, with 
7 = (v3 — l)/4. We find that while N m i n (m) = to for 
< to < 5, but N min (6) > 6. 

We note that the optimal sequences {xj} have been 
defined in terms of H and J, but the polynomial equa- 
tions, and therefore the roots, actually do not depend on 
them at all! A more direct characterization of these has 
not been possible so far. Similar sequences have been 
studied earlier in the context of finding efficient numer- 
ical integration techniques in classical mechanics using 
symplectic integration llj. In quantum Monte Carlo 
studies, the hamiltonian H of the system can often be ex- 



pressed as a sum of two mutually non-commuting terms 
: H = A + B, such that there are efficient procedures to 
apply the operators e* A and e* B to a given state. Then 
using a generalization of the well-known Trotter formula, 
the propagator of the full hamiltonian is approximated 
as 

e tA+tB _ £ xitA e X2tB ^xztA ^XNtB _|_ Q(f ms j (19) 

To make the discretization error small, one wants to 
choose the parameters x\,X2 ■ ■ ■ to maximize to [l2^ . If 
we write — iH = A and — iJHJ = B, then Je~ ltXjH 3 = 
e tXjB , and W N (t) is of the form Ea itT§|) . In the quantum 
Monte Carlo case, it has been proved that for to > 4, one 
cannot find a solution with all {xj } non- negative [13j . In 
our problem, the inversion pulse provides some crucial 
changes of sign in the algebraic equations, making the 
construction of pulse sequences using non- negative {xjYs 
possible for arbitrary to. 

To summarize, in this paper, we have presented an 
algorithm for preserving quantum states using inverting 
pulses, with unequal time interval betwen pulses, which 
does not assume any symmetry properties of the Hamil- 
tonian. The significant improvement over the quantum 
Zcno effect leakage probability seems promising for appli- 
cations to efficient preservation of quantum states against 
decoherence, in quantum computations and communica- 
tions. On the theoretical side, the asymptotic behavior of 
N m i n (m) for large to is an interesting problem for further 
study. 
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